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Abstract. I summarize some constraints on the physics of neutron stars arising from X-ray 
observations of the surfaces of neutron stars, focusing on using models of low-magnetic- field 
neutron star atmospheres to interpret their X-ray spectra. I discuss observations of spectral 
lines, pulsation profiles, X-ray bursts, radius measurements of transiently accreting neutron 
stars in quiescence, crust and core cooling measurements of transiently accreting neutron stars, 
and cooling of young neutron stars. These observations have constrained the neutron star 
mass and radius (and thus the internal composition, and dense matter equation of state), the 
superfluidity and neutrino emissivity properties of the core, and the composition and superfluid 
state of the crust. 



1. Introduction 

The behavior of matter at high densities (several times nuclear density) and relatively low 
temperatures (kT « mc 2 ) cannot be probed in terrestrial experiments. The nonlinear 
aspects of quantum chromodynamics ensure that the equation of state cannot be extrapolated 
unambiguously from well-understood lower-density matter. The interiors of neutron stars (NSs) 
are composed of high-density, relatively cold matter, and the physics of these interiors affect 
the observable properties (e.g. mass, radius, surface temperature) of the NSs. Therefore, 
astrophysical observations of NSs permit constraints on the physics of dense matter, as have 
been summarized in detailed reviews such as [HE]; see also Lattimer (in these proceedings) for 
a complementary discussion. 

NSs can be observed in a variety of manifestations [3], many of which offer possibilities for 
constraints on their internal physics. In this short review, I focus on a few recent advances 
involving thermal X-ray radiation from the surface of the NS [H E]. Other constraints arise 
from, for instance, timing of radio pulses [6l [7] which can provide measurements of the spin rate 
[8] of NSs; neutrino emission from the formation of NSs in supernovae [llj; gravitational waves 
as two NSs spiral together and coalesce [12]; high-energy radiation produced by the decay of 
extremely high magnetic fields in the enigmatic NSs known as magnetars [T3] : and studies of 
binary companion stars and the accretion disks [14J. A key result is NS mass measurements in 
binaries with a radio pulsar, often using tests from general relativity (GR), which suggest that 
the majority of NSs are around 1.4 solar masses (M ) [10J, with some ranging as high as 1.97 

Here I discuss NSs where we observe thermal radiation from the surface (specifically, those 
with low magnetic fields, e.g. B < 10 11 G). As they cool, isolated NSs radiate heat from their 



surfaces which may have nonuniform temperature distributions, leading to flux variations with 
the NS rotation phase. NSs accreting material from a companion star experience occasional 
thermonuclear fusion events on their surfaces, "X-ray bursts", which often provide sufficient 
radiation to briefly expand the NS's atmosphere. Accreting NSs often experience periods of 
quiescence during which material piles up in an accretion disk rather than reaching the NS, thus 
allowing one to study NS surface emission. 

1.1. Thermal surface emission 

NSs are born extremely hot, and their extreme gravity and magnetic fields can accelerate 
particles that heat their surfaces as NSs become older. Thus, many NSs have temperatures 
of order a million degrees over part or all of their surface, and thus glow in X-rays. This X- 
ray radiation is observed by X-ray telescopes on orbiting satellites (as X-rays cannot penetrate 
Earth's atmosphere); the findings reported below have principally come from NASA's RXTE 
and Chandra, and ESA's XMM-Newton satellites. 

The surface radiation from NSs must pass through the thin (few cm) NS atmospheres 
(if present), which modify the radiation through free- free absorption, spectral lines, and 
photoionization jumps. Strong magnetic fields heavily distort the structure of electronic orbitals, 
and thus alter the ionization energies and opacities [15]. Thus, it is easier to get constraints on 
NS physics when studying those NSs which are expected to have relatively low surface magnetic 
fields (£><10 6 T), as the number of free parameters affecting the spectrum is reduced. Accretion 
of material from a companion appears to reduce NS magnetic fields by diamagnetic screening, 
from typical birth fields of 10 8 T to -10 4 T [Tg]. 

Due to the strong surface gravity, NS atmospheres stratify by element within < 100 s fT7] . 
Thus, at times of little or no accretion, the part of NS atmospheres that produce the spectra is 
usually assumed to be pure hydrogen (during X-ray bursts, the composition should be complex, 
and one young NS appears to have a carbon atmosphere [18J). Low-B-fleld hydrogen atmosphere 
models give observable X-ray spectra similar to blackbodies, but shifted to higher energies due to 
the inverse dependence of the electron free- free opacity on photon energy jT9l [20] . The intense 
gravity of NSs shifts the observed spectrum to longer wavelengths, decreasing the observed 
temperature. Although they complicate the spectral interpretation, gravitational effects permit 
constraints on the compactness (mass divided by radius) of the NS. 

2. Constraints 

2.1. Spectral lines 

An unambiguous identification of a clear spectral line due to an atomic transition from a NS 
surface would permit, by measuring how the energy of the line is decreased by escaping from 
the gravitational potential of the NS, a straightforward measurement of the compactness of the 
NS. A tentative identification of narrow absorption lines in a high-resolution X-ray spectrum of 
X-ray bursts from one accreting NS (EXO 0748-676) suggested a gravitational redshift of z=0.35 
[21], consistent with theoretical expectations pQ, and sparking many attempts to find similar 
lines in other accreting NSs. Unfortunately, deeper observations of this NS with multiple X-ray 
telescopes have failed to confirm these lines [22] . Furthermore, the NS rotation period, later 
measured at 552 turns per second [23] , would broaden any lines produced on the NS surface 
[2U[25], implying that the observed absorption lines do not arise from the NS surface. Searches 
for absorption lines in high-resolution spectra of other accreting NSs have also met with no 
success. 

Other spectral absorption lines have been observed in isolated NSs, both in nearby NSs [26|l27] 
and in very young NSs [28j[29]. These spectral lines appear to all be consistent with variants of 
proton or electron cyclotron absorption (or emission) lines [30l|31], thus their study can tell us the 
surface magnetic field strength but not the surface redshift of the NS. Recently, a time- variable 



absorption line has been claimed in a normal pulsar, which shows a magnetic field strength 
incompatible with the line's production by cyclotron absorption [32J. It is difficult to understand 
how a line produced by an atomic transition would have a variable energy. However, the apparent 
lines could instead be due to an intrinsically complicated multiple-blackbody spectrum (cf. [33j), 
which when fit by a single blackbody gives the appearance of absorption lines [32] . 




Figure 1. X-ray pulsation profile of the radio 
pulsar PSR J0030+0451 folded on the NS pulse 
period (repeated twice for clarity), as observed 
by the XMM-Newton telescope and modeled by 
[37] . The two peaks indicate the maximum 
visibilities of the two polar caps. 



Figure 2. Constraints on the mass and radius 
of the NS in PSR J0030+0451 from its pulsation 
profile [37]. The cyan region at upper left is not 
allowed due to causality (the requirement that 
pressure waves must propagate at speeds below 
c). The region ruled out by fits to this NS is 
the green shaded region (at 99.9% confidence); 
the dark blue region is also included at 95% 
confidence. Several loci of mass and radius for 
different NS models are also plotted [38]. The 
dashed line is the canonical NS mass, 1.4 M . 



2.2. Pulsation profiles 

Polar caps of NSs may be heated either by accretion, funneled onto the caps by the NS's 
magnetic field, or by radio pulsar activity, which produces energetic (relativistic) positrons and 
electrons that collide with and heat the polar caps [34J . Either situation produces temperature 
inhomogeneities over the NS surface, which combine with the angular dependence of intensity 
from the atmosphere and ray-tracing in a general relativity metric to produce a complicated 
pulse profile [35J. Generally, the Schwarzschild metric combined with appropriately oblate (due 
to their rotation) neutron stars, Doppler effects, and light-travel time delays provides a good 
approximation to the full GR lightcurve calculation [36] . 

Recent work has focused on accreting millisecond pulsars and radio millisecond pulsars, which 
are both believed to have relatively low magnetic fields. The accreting millisecond pulsars, 
due to continuing accretion of a range of elements, may have spectra reasonably described by 
blackbodies, while radio millisecond pulsars, without continuing accretion, should typically have 



hydrogen atmospheres. The simplest constraint from modeling the pulsation profiles is that 
a relatively sharp, peaked profile requires a larger star, since radiation from a more compact 
star suffers more gravitational bending and thus produces a smoother pulse profile. A strong 
constraint on the radius of the NS in the radio pulsar PSR J0030+0451 was calculated by [37] , 
giving R > 10.4 km (99.9% confidence) for an assumed mass of 1.4 M , due to its highly peaked 
pulsation profile (Fig. 1). 

More complicated modeling of the emission from accreting millisecond pulsars (which have 
components of their spectra that do not come from the surface, e.g. Comptonized radiation 
from above the surface, and reflected light) produces constraints in mass and radius space. For 
the well-studied object SAX J1808. 4-3658, [39J find a relatively low mass below 1.5 M (using 
a reasonable constraint on the inclination). Although this modeling involves more parameters, 
the mass constraint is consistent with optical studies of the companion star [40J. 
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Figure 3. Flux vs. temperature for X-ray 

spectra taken at various times during multiple X- 
ray bursts by the same NS [44J. Diagonal lines 
indicate the best-fit blackbody normalization and its 
uncertainty. The majority of bursts are consistent with 
one normalization. 
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Figure 4. Inferred blackbody nor- 
malization vs. X-ray flux for the X- 
ray burst spectra of Fig. 3 [44J . 



2.3. X-ray bursts 

Accreting NSs often experience unstable nuclear burning on their surfaces, leading to rapid X-ray 
brightenings known as "bursts" [H] . Several authors have calculated models for NS atmospheres 
during bursts, showing how the measured temperature (defined by fits to a blackbody spectrum) 
relates to the effective surface temperature; this relation is parametrized as a "color correction 
factor" [151 S6]. The burst may expand the NS atmosphere; the radiative flux necessary to 
overcome gravity and push the atmosphere away, the "Eddington limit", depends on the NS 
mass, but the observed flux will be decreased by the NS's gravity, and thus its compactness. As 
the radiation from the burst comes from the NS surface (rather than the physically complicated 
accretion flow that affects the radiation at other times), measurement of the temperature and 
flux of these bursts (and knowledge of their distances) determines the emitting area, and thus 
(after including GR corrections) the NS radius [42]. Similarly, measuring the radiation flux 
during bursts that expand the atmosphere tells us about the surface gravity of the NS [43] . 

Using time-resolved spectra of several bursts showing atmospheric expansion, [47l [48j [49] 
calculated constraints on the mass and radius of three NSs, finding relatively small NSs 



(with radii of 8-11 km). These calculations have been criticized by [50] on their choices 
of color correction factors, distances, and surface chemical compositions, and for not being 
internally self-consistent; the key problem is that the best-fit values of the observable properties 
produce complex (non-real-valued) values for the NS mass and radius. [50] present alternative 
calculations tending toward slightly larger (^12 km) NS radii, favoring the interpretation that 
the observed radius of the atmosphere remains larger than the true NS radius throughout the 
burst. [51] use their color correction factors to derive a large radius (>14 km) for a burst from 
one NS, though they have also been criticized [2] for the relatively poor fits of their spectra to 
blackbodies. 

Several uncertainties surround the interpretation of X-ray bursts. A critical question is 
whether the NS radii measured during the course of a burst, and between bursts, are identical 
and equal to the NS radius [52l [44] . Fig. 3 (from [44]) shows temperature vs. flux at different 
times across multiple bursts from one NS, showing consistency of most burst spectra with one 
inferred radius (Fig. 4), which suggests that this uncertainty is manageable. Other questions 
include whether the burst fluxes consistently reach the same Eddington limit [53} IM] , whether 
they show significant anisotropy [55] , and whether they vary in their chemical composition [56] . 
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Figure 5. Constraints on the 
mass and radius of two transiently 
accreting NSs at known distances, 
as inferred from fitting their X-ray 
spectra with hydrogen atmosphere 
models. Hatched areas indicate the 
90% confidence regions for NSs in 
47 Tuc (red) [58] and M13 (blue) 
[59] , with outside contours (solid for 
47 Tuc, dashed for M13) marking 
the 99% confidence regions of each. 
Some sample NS equations of state 
are indicated by dotted black lines 
[38] . The 90% confidence contour 
regions for these two NSs do not 
overlap. 



2.4- Quiescent accreting neutron stars: radius 

Many accreting NSs show transient behavior, with periods of little or no accretion, in which the 
surface of the heated NS may be studied. Assuming a low-magnetic-fleld hydrogen atmosphere 
without anisotropics, and a distance known by other means (e.g. by studying the other stars 
in its stellar cluster), tight constraints may be placed on the NS mass and radius by measuring 
the temperature and flux, and thus inferring the radius [57] . Due to the gravitational redshift, 
curved constraints in mass and radius are produced (e.g. Fig. 5). 

Deep X-ray observations of globular clusters containing transient accreting NSs have placed 
tight constraints on the mass and radius of several NSs [581 E3 EQl [61] . Not all these constraints 
are consistent with each other (see Fig. 5), which may be due to instrumental effects [58] , and/or 



variations in atmospheric composition between NSs [61J. Combining results from these NSs with 
X-ray burst data (see above) gives robust constraints on the NS radius (between 10.4 and 12.9 
km for 1.4 M NSs), and thus on the dense matter equation of state [62J. 
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Figure 6. Cooling curves for 
various neutrino emission scenarios 
[65], compared with measurements 
or upper limits of the NS luminos- 
ity outside accretion, vs. the time- 
averaged mass transfer rate [67J. 
The dotted line relates heat input 
(via mass transfer) to the NS lumi- 
nosity between accretion episodes 
for "standard" (slow) cooling, while 
the other diagonal lines suggest 
the maximal cooling achievable by 
specific enhanced cooling processes. 
The effect of a factor of 1.5 error in 
distance is indicated. SAX J1808- 
36's cold temperature (D) rules out 
most neutrino emission scenarios. 



2.5. Quiescent accreting neutron stars: core temperature 

Observations of transient NSs between periods of accretion measure the temperature that the 
NS has been heated to during accretion episodes. The kinetic energy of the infalling material 
is radiated away immediately, but the added material compresses the deep crust, leading to 
nuclear reactions that deposit heat deep in the crust, some of which flows into the core [631164]. 
Averaged over long periods including many accretion episodes, the heat flowing into and out 
of the core will reach equilibrium, and the NS surface (between outbursts) will reflect the 
core temperature. Neutrino emission from the NS core dominates the NS's energy loss, either 
through "standard" cooling (neutron-neutron neutrino bremsstrahlung processes) or enhanced 
neutrino cooling mechanisms such as direct URCA (involving protons or hyperons) or URCA- 
like processes involving (e.g.) pion or kaon-like Bose-Einstein condensations in the core, which 
will decrease the NS surface temperature J65J [66] • Differences in the surface temperature of 
NSs with the same average mass transfer rate (and thus heating history) can be explained by 
different NS masses, as this will change the central NS density and perhaps composition, and 
thus alter which cooling mechanisms may be accessed. 

Strong evidence for enhanced neutrino cooling of the cores of transient NSs was provided 
by relatively cold NSs whose accretion histories suggested higher NS temperatures [68} \69\ 170]. 
Perhaps the clearest evidence for enhanced neutrino cooling is the cold NS SAX J1808. 4-3658, 
where the mass transfer rate is very well-known (Fig. 6); this object may require direct URCA 
processes involving protons or hyperons [7T|l67]. An important question is whether the observed 
accretion history accurately reflects the long-term averaged mass transfer rates (see [72] and [73] 
for differing views on this NS). 
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Surface temperature 
for the transiently 
accreting NS MXB 1659-29 after 
accretion stopped, vs. models 
with different choices of crustal 
impurities (the initial temperature 
after accretion is adjusted to match 
the first point) [80]. X-axis is 
days since accretion stopped, y-axis 
is NS surface temperature. The 
level of impurities (elements that 
disrupt the crystalline structure) in 
the crust is tightly constrained. 



2.6. Quiescent accreting neutron stars: crust temperature 

Observations of transient NSs immediately after long periods of accretion allow measurements 
of the cooling behavior of the crust, which radiates the heat produced by nuclear fusion within 
the crust within a few years, coming back to equilibrium with the core [74J. Crustal cooling has 
now been tracked for five transients over several years [T5J, [T6J, [771 EH1 [79] . 

The crustal cooling behavior is dominated by the crust's thermal conductivity, which is set 
by the level of defects in the crystalline crust, and its heat capacity, which is strongly reduced 
if the neutrons are superfluid. Thus, measuring crustal cooling constrains the composition of 
the crust, and thus how nuclear reactions occur in the crust [81]. Comparisons of the crust 
cooling curves with models (Fig. 7) show that the crust is thin, highly conductive (so crystalline 
with low impurities), and that the neutrons are indeed superfluid |82t [80] . A key remaining 
uncertainty is whether the measured NS surface temperature is affected by continuing, low-level 
accretion [77], which might appear to increase the rate of the temperature decay. 



2.7. Young cooling neutron stars 

The thermal behavior of NSs can also be studied by measuring the temperatures of young NSs of 
known ages, from 300 to roughly a million years old [65J. Here there is less uncertainty about the 
heat input and timescales (compared to transiently accreting NSs), but there is also a smaller 
range of temperatures. Allowing for different compositions of the outer layers of the crust, 
[83] showed that the known young NSs were consistent with a "minimal" (standard) cooling 
model, without the enhanced cooling mechanisms suggested to explain some of the transient 
NSs. However, selection effects may affect this, as identifying young NSs often requires relatively 
high temperatures. There are several supernova remnants which may contain young NSs, but 
sensitive, unsuccessful X-ray searches for them have constrained any such NSs to relatively low 
temperatures. Considering the young ages of these supernova remnants, such cold NSs so early 



Figure 8. Cooling curves 

(solid lines) for NSs with a range 
of masses (varying other param- 
eters can reach the blue dashed 
lines), compared with observations 
(crosses) of the temperatures of 
young NSs of different ages. The 
drop in model curves at around 100 
years is due to the neutron super- 
fluid transition. The slope of the 
model curve passing through the 
Cas A NS data point (red) also 
matches the observed temperature 
drop of the Cas A NS [85]. 

may require enhanced neutrino cooling [83] . 

The youngest known NS, in the 330- year-old supernova remnant Cassiopeia A (Cas A), shows 
evidence of rapid cooling, by ^3% over 10 years [86]. It is too hot now to have experienced 
enhanced cooling, as those mechanisms would operate throughout its history and make it very 
cold already. But its current temperature drop is too fast for normal cooling, requiring a 
recent transition to more rapid cooling. This can be explained if the neutrons in the NS core 
are currently transitioning to a superfluid state; as the neutrons bind into the quantum pairs 
required to enter the superfluid state, they release the "gap" energy of the pairing in the form 
of neutrinos [87]. Models of the temperature decline from such a transition match the data, for 
a critical temperature of the neutron superfluid of 5 — 9 x 10 8 K [88j [85]; see Fig. 8. The key 
question here is how well the calibration of the X-ray telescope (always a complex question) is 
known. This is the first observationally driven estimate of the superfluid critical temperature 
for any superfluid other than the helium-3 and -4 isotopes on Earth, and a beautiful illustration 
of the potential NSs can offer for constraining properties of the superdense matter. 
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